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Abstract 

Nonextensive aspects of the degree distribution in Watts-Strogatz (WS) small-world net- 
works, Psw{k), have been discussed in terms of a generalized Gaussian (referred to as 
Q -Gaussian) which is derived by the three approaches: the maximum-entropy method 
(MEM), stochastic differential equation (SDE), and hidden- variable distribution (HVD). 
In MEM, the degree distribution Pqik) in complex networks has been obtained from 
Q-Gaussian by maximizing the nonextensive information entropy with constraints on av- 
erages of k and k"^ in addition to the normalization condition. In SDE, Q-Gaussian is 
derived from Langevin equations subject to additive and multiplicative noises. In HVD, 
Q-Gaussian is made by a superposition of Gaussians for random networks with fluctuat- 
ing variances, in analogy to superstatistics. Interestingly, a single Pqik) may describe, 
with an accuracy of | Psw{f^) ~ Pq{^) I~ 10~^, main parts of degree distributions of 
SW networks, within which about 96-99 percents of all k states are included. It has 
been demonstrated that the overall behavior of Psw{k) including its tails may be well 
accounted for if the fc-dependence is incorporated into the entropic index in MEM, which 
is realized in microscopic Langevin equations with generalized multiplicative noises. 
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1 INTRODUCTION 



In the last several years, there is an increased interest in two subjects of statistical physics: 
(1) nonextensive statistical mechanics [l]-[4] and (2) complex networks [5]- [10]. As for the 
first subject (1), Tsallis has developed the nonextensive statistics (NES), proposing the 
generalized entropy given by [1][2] 



where q denotes the entropic index and pi the probability distribution of the state i. It 
is noted that in the limit oi q — 1, Sq reduces to the entropy in the Boltzman-Gibbs 
statistics given by 



The NES has been successfully apphed to a wide range of nonextensive systems including 
physics, chemistry, mathematics, astronomy, geophysics, biology, medicine, economics, 
engineering, linguistics, and others [4]. 

The subject (2) has been initiated by the seminal paper of Watts and Strogatz [5] . 
Since then considerable study has been made on complex networks [9] [10]. Existing and 
proposed networks are classified into two types: 

(i) Networks where the number of nodes is constant and where the degree distribution 
P{k) for a given node to be connected to k neighbors, has a maximum at or near < k > 
with finite width where < ■ > denotes the average over P{k). This type includes random 
[11] and small-world networks [5] [6]. 

(ii) Networks where the number of nodes is not stationary and where the distribution varies 
between the exponential and power forms. This type includes scale-free (SF) networks, 
which were originally proposed by Barabasi and Albert [8] with a growth of nodes and their 
preferential attachments. Many models and mechanisms have been proposed not only for 
growing but also for non-growing SF networks with geographical and non-geographical 
structures. 

It is possible that the two subjects (1) and (2) are intimately related [12]-[16]. By 
using the nonextensive information entropy as given by Eq. (1), Wilk and Wlodarczyk 
have shown that degree distribution P{k) of SF networks which belongs to type (ii), is 
given by [12] [13] 



(2) 




(3) 
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where fi =< k > and < ■ > denotes the average over P{k). The expression (3) is 
universal because in the hmit of g = 1 it reduces to the exponential distribution given by 
P{k) oc e"*^/'^, and for k ^ fJ'/iQ — 1) it yields the power law given by P{k) oc k^'^^'^^'^^\ 
Thus the expression (3) successfully accounts for degree distribution in SF networks. 
Scares et al. [14] have proposed the two-dimensional growth model for SF networks 
where the degree distribution is given by Eq. (3) with and q expressed in terms of 
model parameters. Thurner and TsaUis [15] have discussed scale-free gas-like networks, 
showing their nonextensive aspects. It has been pointed out by Abe and Thurner [16] 
that degree distribution in scale-free networks may be generated by superimposing Poisson 
distributions relevant to random networks [11] with a proper hidden- variable distribution, 
just as in superstatistics [17]. 

As for the type (i) networks, degree distribution P{k) of random network is given by 
the binomial distribution [11]: 

PB{k) = cf- V (1 - c)^-l-^ (4) 

where c = ji/N and = m\ / n\ {m — n)\. In the limit of large iV, Eq. (4) becomes the 
Poisson distribution given by 

k 

Pp{k) ^ Pp{k; f,) = |-e-^ (5) 

On the contrary, P{k) of Watts-Strogatz (WS) small-world (SW) networks [5] is given by 
[18] 

n=0 

(for k>n/2 and large N) (6) 

where f{k,fi) = min(A; — /i/2,/i/2) and p (= 2Nr/Nfi) denotes a randomness parameter 
for random rewirings of Nr edges from a regular network [5]. Figure 1 shows Pswik) with 
N = 100 and fi = 10 for various randomness p [19]. Results shown in Fig. 1 are obtained 
by simulations with 1000 trials after Ref. [5], because the analytic expression given by Eq. 
(6) is valid only for k > n/2 and large N. For p = (regular networks), Psw{k) becomes 
a delta-function at k — /i. With more increasing p, Psw{k) has a peak at k — jj, with 
a wider width. As a comparison, we show by open circles in Fig. 1, Psik) of random 
networks given by Eq. (4). Note that even for p — 1, Psw{k) of the WS-SW network is 
different from Pb(^) of random networks because the former is not locally equivalent to 
the latter with some memory of the starting regular network [18]. 
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We suppose that Psw{k) shown in Fig. 1 may be described or approximated by a 
simple function Uke Eq. (3) if the nonextensive property of SW networks is properly 
taken into account. The purpose of present paper is to investigate such a possibility. This 
is indeed possible in an approximate sense, as will be demonstrated in the present paper. 

The paper is organized as follows. By using the maximum-entropy method (MEM) 
with constraints on < A; >= /i and < >— p besides the normalization condition, we get 
in Sec. 2, the generalized Gaussian with a maximum at k — /i (hereafter referred to as 
Q-Gaussian). In contrast, a conventional ^-Gaussian distribution which has been widely 
employed in many subjects such as self-gravitating stellar systems [20] and ideal gas [21], 
has a peak at the zero value. In Sec. 3, we discuss the stochastic-differential equation 
(SDE) by using Langevin equation with additive and multiplicative noises. In Sec. 4, SW 
networks are constructed by hidden- variable distribution (HVD). An analysis of degree 
distributions of SW networks is made in terms of Q-Gaussian in Sec. 5. Section 6 is 
devoted to more detailed analysis, including the /c-dependence in the entropic index in 
MEM, which is shown to be realized in Langevin equations with generalized multiplicative 
noises. Our conclusions are presented in Sec. 7. 

Before going to Sec. 2, we summarize various distributions used in this paper as 
follows. PsfO^), SF networks; PswO^), SW networks; PBik), binomial; Pp{k), Poisson; 
Paik), Gaussian; p^, Q-Gaussian; Pgik), an escort probability of p^. 

2 Maximum-entropy method 

2.1 Shannon information entropy 

We have adopted a network with N nodes whose adjacent matrix is given by — Cji — 1 
for a coupled pair (i, j) and zero otherwise. The degree distribution P{k) for a given node 
to have k neighbors is defined by 

^(^) = ^E<'^(^-E%)>G, (7) 

» 3 

where <>g denotes the average over graphs. The averaged coordination number // is 
given by 

= ^EE<Ci^>G, (8) 

* 3 

= E^(^)^- (9) 

k 
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The coupling connectivity R is defined by [19] 

^ = ^EEE<Ci.-Ci^>G, (10) 

= -,T.P{k)e. (11) 
k 

The couphng connectivity expresses a factor for a cluster where the two sites j and 
£ are coupled to the third site i, but the sites j and £ are not necessarily coupled. In 
contrast, the clustering coefficient C defined by 

^ = TTT^EEE < >G, (12) 

^^/^ i J e 

expresses a factor forming a cluster where the three sites i, j and £ are mutually coupled 
[22]. Unfortunately C cannot be expressed in terms of P{k). 
First we consider the Shannon entropy given by 

S ^ - Y,Pik)^^P{k). (13) 

i 

Assuming that and R of the adopted networks are given, we obtain P{k) by MEM with 
the three constraints given by 

= 1, (14) 

k 

j:p(k)k = (15) 

k 

Y,P{k)k^ = fi^R = p\ (16) 

k 

We get 

P{k) oc e-^^''+^'''\ (17) 

where (3 and 7 are Lagrange multipliers relevant to the constraints given by Eqs. (15) 
and (16), respectively. 

2.2 Nonextensive information entropy 

Next we extend our discussion by employing the generalized entropy first proposed by 
Tsallis [1, 2]: 

S, = ^f^, (18) 
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which reduces to the Shannon entropy in the hmit of g = 1. Here pk denotes the probabihty 
distribution for couphngs of k neighbors whose exphcit form will be given shortly [Eq. 
(24)]. The three constraints corresponding to Eqs. (14)-(16) are given by 



k 

k 



1, 



with 



Pk 



P{k) = ^ 



k 



(19) 
(20) 
(21) 

(22) 
(23) 



where P{k) is the escort probability expressing the degree distribution [2]. By maximizing 
the entropy with the three constraints, we get 

'1-g^ 



Pk oc 



1 - 



1-9 



Equation (24) may be rewritten as 



Pk = ^expg 



with 



Zq = 


J2 exp 




k 


C = 


m 












27' 



1 + 



1 -1 



(24) 
(25) 

(26) 

(27) 
(28) 



where expg(a;) expresses the g-exponential function defined by 

exp^(a;) = [1 + (1 - g)a;]^, for 1 + (1 - g)x > (29) 

= 0. for l + {l-q)x<0 (30) 

Parameters of ( and rj, newly introduced in place of /3 and 7, are determined by the 
constraints given by Eqs. (20) and (21) as 

= j:p(k)k, (31) 

k 

= T.Pik)k', (32) 

k 
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with 



(expJ-C(A;-77)2])« 



(33) 



Efe (exp, [-C{k-r]r]y 
which are numerically evaluated in general. 

2.3 Gaussian 

In order to get some analytical results, we hereafter assume that /c is a continuous variable 
varying from — oo to oo, which may be justified for a. Equations (31)-(33) yield 



2 2 
p - fi , 



:i-5)c 



1 

2C' 



^\2-> g-l 2/ 



for g < 1 
for g = 1 
for g > 1 



L(?-i)Cj 

where B{x,y) denotes the beta function. Equations (35)-(38) lead to 

1 



c 



with 



Prom Eqs. (25), (34) and (39), we get 



3-g 



1 



Pk=Pk 



with 



D{q,r) = 



27raD{q,l) 



2r 



exp„ 



l_q \2r + l-q) lr(^ + i)r 



for g < 1 
for q = I 
for g > 1 



(34) 

(35) 

(36) 
(37) 
(38) 

(39) 
(40) 

(41) 

(42) 
(43) 
(44) 



V r^^) 

r{x) being the gamma function [23]. Prom Eqs. (22), (23) and (41), we get degree 
distribution given by 

P{k) = Pgik) = ^ ^ I exp,, 



27r a D(q,q) 



{k - ixf 



2i/(72 



(45) 



p2 in Eq. (41) expresses the generalized Gaussian which is called Q- Gaussian in this 
paper. It is easy to see that D{q,r) reduces to unity in the limit of g = r = 1 because 
lim|2|^oor(^ + a)/r{z)z°' = 1 [24], then both and Pqik) reduce to Gaussian: 

Poik) = PG{k;fi,a) = e-^, (46) 



with the average fi and variance a^. Equations (39) and (40) show that q < 3 because 
> for the positive definiteness of fluctuations. We note from Eqs. (29) and (30) that 

2q 

Pqik) oc {k — fj) 1-^ , for \ k — fj, \ / a ^ rc, q > 1 

= 0, for I - /X I /cr > Tc, g < 1 (47) 

with 



/ 3- 


q 


1 Q- 


1 \' 



(48) 

It is worthwhile to note the difference between Q-Gaussian and the conventional q- 
Gaussian which is given by p(-u) oc exp^[— derived from Eq. (3) with a replacement of 
k (ii^/2) [20] [21], u standing for the generahzed velocity and the relevant Lagrange 
multiplier. Alternatively, g-Gaussian is obtained from Q-Gaussian with = in Eq. (41). 

2.4 Model calculations 

In order to first get a broad insight to the Q-Gaussian, we show in Figs. 2(a) and 2(b), 
some numerical results of and Pqik) calculated by using Eqs. (41)- (45). Figure 2(a) 
shows with = 10 and o" = 1 for various q values of g = 0.5, 1.0, 1.5, 2.0 and 2.5. 
The Q- Gaussian is fat-tailed for g > 1 and compact support for g < 1. Fig. 2(b) 
shows Pqik) with yU = 10 and a = 1 for various g values corresponding to Fig. 2(a). We 
note that a profile of Pqik) is rather different from that of p^: both are normalized by 
/ ci/cpfc = / dk Pqik) = 1. For g = 1.0, and Pq{k) reduce to the Gaussian distribution 
located at k — /i — 10. With more increasing g above unity, Pq{k) has a sharper peak 
with narrower width. In contrast, when g is more decreased below unity, Pq{k) has a 
wider width but has no magnitudes at | A; — // |> Tc cr, for example, | k — /i \> 2.236 for 
g = 0.5. 
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3 Stochastic differential equation approach 
3.1 Formulation 

Recently Anteneodo and Tsallis [25] have studied the Langevin equation subject to addi- 
tive and multiphcative noises, which lead to nonextensive distributions. We here assume 
that Langevin equations are given by 



dk- 

= -A {k + g{h) Ut) + rii{t), {i = l-N) 



(49) 



where ki{t) denotes a real variable expressing the number of couplings of the node i: A 
the relaxation rate: g{k) an arbitrary function whose explicit form will be given later 
[Eqs. (52), (88) and (89)]: Tjiit) and {j(t) are additive and multiplicative white noises, 
respectively, as given by 



< iiit) i^it') > = 2M6ij6{t-t'), 



(50) 
(51) 



means of noises being vanishing. Equation (49) shows that the number of couplings ki is 
influenced by both additive and multiplicative noises. In particular, effects of multiplica- 
tive noises given by g{ki) depend on the present state of node ki{t). 
When g{k) is given by 

g{k) = k-fi, (52) 
we get the Fokker-Planck equation for p{ki, t) (= p{k, t)) given by [25] 



(53) 



Its stationary solution is given by [25] [26] 



M- 



Pk =p(A;,oo) oc 
Equation (54) may be rewritten in two ways as 



A+M 
2M 



(54) 



Pk oc 



oc 



{k-J^f 

2a' 



2l/(72 



1-9 



l-q 



(55) 
(56) 
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with 

A + 3M 

= vAt' (58) 
A + M ^ ^ 

= p (59) 

where is given by Eq. (40). Note the difference between a in Eq. (58) and a in Eq. 
(59). The Q-Gaussian given by Eq. (55) with ^ = is adopted in [25]. Equation (56) 
agrees with the result obtained by MEM given by Eq. (41) for q > 1. 

Equation (57) shows that the nonextensivity arises from muhiphcative noises. In 
the case of no multiplicative noises (M = 0), we get q = 1, leading to the Gaussian 
distribution: 

Pk = Paik) = e-^. (60) 

v27r(T 



3.2 Model calculations 

We show model calculations, simulating Langevin equation given by Eq. (49) with N — 
100: the Heun method is employed with a time step of 0.0001. Figures 3(a) and 3(b) 
show time courses of ki{t) [= k{t)] for {A,M) = (1.0,0.0) and (1.0, 0.5), respectively, 
with g{k) = k — fi, X = 1.0 and /i = 10. The number of k, starting from an initial value of 
/i = 10, fluctuates around k = fihy effects of noises. A comparison between Figs. 3(a) and 
3(b) shows that fluctuations are much increased by effects of multiplicative noises. Circles 
in Figs. 4(a) and 4(b) show degree distributions P{k) for (A, M) = (1.0,0.0) and (1.0, 
0.5), respectively, obtained by averages over N — 100 for 50 < t < 100 [27]. Solid curves 
in Figs. 4(a) and 4(b) express Pgik) given by Eq. (45) for (g, a) = (1.0, 1.008) and (1.667, 
0.984), respectively: adopted q values are calculated by Eq. (57) and a are evaluated by 
using Pq(A;) obtained by simulations [27]. Degree distribution for M = in Figs. 4(a) 
shows Gaussian. When the multiplicative noises with M — 0.5 are introduced, degree 
distribution P{k) has tails expressed by the power law as shown in Fig. 4(b). Results of 
stochastic differential equation are in good agreement with those of Q-Gaussian shown by 
solid curves. 
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4 Hidden-parameter distribution approach 



Employing Poisson distribution Pp{k) for random networks, Abe and Thurner [16] pro- 
posed the hidden-parameter expansion to get 

/ dpi' Pp{k;pt')U,{pL')^P{k), (61) 
Jo 

where ni(//) stands for distribution of the hidden variable /i. For scale- free distribution 
PsF{k) given by 



they obtained [16] 



with 



PsF{k) = A{k + ko)-\ (62) 



ni(/i) oc ji ^ , for large ji (63) 



5 = 7, (64) 



where 7 > 1, and A and ko denote constants. 

It should be noted that degree distribution of random networks may be described not 
only by Poisson Pp{k) but also by Gaussian distribution Pcik). Pp{k) is obtained for 
large N from binomial distribution Psik) as shown by Eq. (5). Alternatively, Pcik) may 
be derived from binomial distribution when k in Ppik) is assumed to be a real variable 
as follows. Expanding Psik) around k = /i = Nc where Psik) takes the maximum value, 
we get 



1 rfc-,i)2 

^ - Poik), (65) 



. 27r(7 

where = Nc{l — c). 

When generalizing the hidden-parameter method proposed by Abe and Thurner [16], 
we may assume that P{k) is expressed as a superposition of Gaussians Paik) for random 
networks, as 

poo roo 

/ d/i' / da' Paik; /i', a') U{ii', a') = P{k), (66) 
Jo Jo 

where n(/x, a) denotes distribution for hidden variables /j, and a. When Pswik) obtained 
by simulations or an analytical expression given by Eq. (6) is substituted to P{k) in Eq. 
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(66), n(/i, o") will be determined by inversely solving Eq. (66). This would be interesting 
but seems rather difficult. 
When we assume 

U{f,',a')=U2{a')d{f,'-fx), (67) 

Eq. (66) becomes 

/ da' PG{k;ii,a')U2{a')^P{k), (68) 
Jo 

where n =< k >. We may get 112(0"') for Psw{k) given by Eq. (6), in principle, from Eq. 
(68), which is, however, problematic because Pswik) is valid only for k > n/2. 
Instead, we have obtained 112 (cr) for Q-Gaussian from Eq. (68) by 

/ da' e Ha a' = exp^ ^ V 

Jo V2^a' ^ ' ^[ 2va^ 

where v — {3 — q)/2 [Eq. (40)], and normahzation factor inp^ is neglected for a simplicity 
of calculation. By using the change of variables: l/2a''^ — a, l/2ua'^ — ckq, and {k — jiY — 
X, we get 



(69) 



with 



/ da e~°'^ f{a) = expg{—aox), (70) 

J 



/(.) = J,f±)nia)"\ (71) 



^ \2aJ \\2a 
1 / n 



r(f) V2ao 



2 n 1 



q;2-1 e 2ao, (72) 



ao = = E{a), (73) 

(74) 

where /(a) in Eq. (72) denotes the F- (or x^-) distribution function of the order n, £"(«) 
the average of the variable a over /(a), and £'(a^) its variance [28] [29]. Equation (70) is 
nothing but the Laplace transformation. By using the relation derived from Eq. (71): 

n.(.) = (^) / (^) . (75) 

we get 

n2((7) oc a-\ for large cr (76) 

with 

^ = 7=^. (T7) 
12 



where 7 denotes the index for the expression given by 

p2 oc (k-fj,)-'^. for large A; (78) 

Eqs. (76)- (78) should be compared to the result of Ref. [16] given by Eqs. (62)- (64), 
related discussion being given in Sec. 6. 

Figure 5 shows Log plots of 112 (u) for Q-Gaussian with various q values against 
c/(co-\/^) where ctq = I/v^Oq for q — 1.0. In the limit of g — > 1.0 {n — > 00), 112(0") 
becomes the delta function at a — ao. As 5 is more increased than unity, 112 (cr) has a 
wider width, expressing larger fluctuations in the variable a. For q — 1.5, 2.0 and 2.5, 
the index 5 given by Eq. (77) is 4, 2 and 1.33, whose slopes are shown by chain curves in 
Fig. 5. 

If we express Gaussian in Eq. (66) as 

PG(k;^i',a')^ij(^^-^y (79) 

we get 

P{k) = dfi' da'^[^^^ n(/.',a'), (80) 

which is similar to the wavelet transformation [30]. It is interesting to obtain analytical 
expression of n(/x, a) for a general P{k). 

5 Analysis of degree distribution in SW networks 

We will discuss an apphcation of the Q-Gaussian to an analysis of Psw{k). The coupling 
connectivity R has been calculated with the use of Eq. (10) by simulations (1000 trials) 
for SW networks with = 100 and /i = 10 as a function of p [19]. We have obtained a 
by using a = fiy/ R — 1 derived from Eqs. (21) and (35). Chain and dashed curves in Fig. 

6 show the p dependence of R and a, respectively. With increasing p, both R and a are 
monotonously increased. 

We have tried to fit Pqik) in Eqs. (41)-(45) to Pswik) obtained by simulations, 
choosing a proper q value with and a calculated for a given p shown in Fig. 6. When 
we have chosen the q value such that maxima of the two distributions agree, their whole 
shapes are in good agreement. Figures 7(a)-7(f) show such fittings for p — 0.1, 0.2, 0.4, 
0.6, 0.8 and 1.0. In the case of p = 0.1 shown in Fig. 7(a), Pq(A;) with q = 1.35, a = 0.97 
and jJL — well reproduces the relevant degree distribution plotted by circles. For a 
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comparison, wc show also results with the use of g = 1.5 and 1.2 whose agreements are 
worse than that of g = 1.35. Gaussian (g = 1.0) for /x = 10 and a = 0.97 is in significant 
disagreement with Psw{k) aX 9 < k < 11. In the case of p = 0.2 shown in Fig. 7(b), 
Pqik) with q = 1.2, a = 1.34 and /x = 10 well explains Psw{k) obtained by simulations, 
results with g = 1.0 and 1.4 being also plotted for a comparison. In cases of p — 0.4, 0.6 
and 0.8 shown in Figs. 7(c), 7(d) and 7(e), we have obtained a good agreement between 
Psw{k) and Pqik) with g = 1.0 expressing the Gaussian distribution. In contrast, in the 
case of p = 1.0 shown in Fig. 7(f), the result of g = 0.95 yields a better fit than that of 
g = 1.0. 

At a glance at Figs. 7(a)-7(f), we have an impression that Psw{k) is well described 
by Pqik). It is, however, not true in a strict sense, because at | A; — // | /u ^ Tc, Psw{k) 
follows the exponential function while Pqik) obeys the power law (g > 1), as Eq. (47) 
shows. This is clearly seen in Fig. 8(a)-8(f) where results in Figs. 7(a)-7(f) are shown in 
the Log plots. For example, in the case of p = 0.1, the agreement of Pqik) with Psw{k) is 
good for I A; — |~ 3 as demonstrated in Fig. 7(a), but it becomes worse for | /c — |~ 3. 
Nevertheless, it should be stressed that main parts of Psw{k) are well described by a 
single Pq{k) as shown in Fig. 7(a)-7(f). When Psw{k) and Pqik) are in good agreement 
at /c G [/Li — fca, yU + ko\ with an accuracy | Pswik) — Pq{k) \< 10~^, the number of k states 
included in this region is given by 



which yields = 0.992 for = 3 in the case of p = 0.1. Similarly, we get (p, ka,na) — 
(0.2,4,0.993), (0.4, 4, 0.975), (0.6, 5, 0.983), (0.8, 5, 0.970), and (1.0, 5, 0.959). 

6 Discussions 

We will present more detailed analysis, trying to reduce the discrepancy between Psw{k) 
and Pq{k) at their tails in this section. Psw{k) calculated by simulations for SW networks 
with p = 0.1, /x = 10 and N — 100 is again shown by filled circles in Fig. 9(a), where 
the analytical result vahd for large N given by Eq. (6) is also shown by filled triangles 
for a comparison. We note a deviation of Psw{k) from Pq{k) shown by solid curve at 
I /c — /X |~ 3, where Psw{k) obeys the exponential (Gaussian) law. This suggests that p^ 
in SW networks behaves as 



/ dk P{k), 



(81) 



Pk oc expg 



(3-g)<72 



for small \ k — jj, 
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oc e 2a'^ . for large \ k — /i] (82) 

If we incorporate the /c-dependence into the entropic index, the entropy in Eq. (18) is 
generahzed as 

V n** — 1 

S, = (83) 

where denotes the /c-dependent entropic index. MEM with the constraints given by 
Eqs. (19)- (21) with a replacement of g — > yields 



with 




from which P{k) is given by Eqs. (22) and (23) with a replacement of g ^ g^. Wc expect 
that the result given by Eq. (82) may be realized in Eq. (84) if g^. = g for small \ k — /j, \ 
and gfc = 1 for large \ k — /j, \. 

Bearing these consideration in mind, we have assumed the /c-dependent entropic index 
given by 

Qk = i + {q-i)e{k- i^ + kuO) e{i^ + k^-k,e), (86) 

with 

= ;PWTT' ^^^^ 

where ki denotes the characteristic k value. Equation (86) shows that g^ = g for | A; — // |~ 
ki and g^ = 1 for | A; — |~ /ci. Open squares in Fig. 9(a) show P{k) calculated by using 
Qk given by Eqs. (86) and (87) with ki — A and 9 — 1. The /c-dependence of P{k) changes 
from Q-Gaussian to Gaussian with increasing | k — /i | , and it is in better agreement with 
Psw{k) than Pq(A;). 

In order that the /c-dependence in the entropic index given by Eq. (86) is realized in 
the Langevin model, we have tried to modify a function g{k) for multiplicative noises. 
Open circles in Fig. 9(b) express the results of simulations of Langevin equation (49) with 
g{k) = k — fi given by Eq. (52) for A = 1, A = 0.941 and M = 0.212, which are chosen 
from Eqs. (57) and (59) with g = 1.35 and a = 0.97 for p = 0.1 [27]. In contrast, open 
squares in Fig. 9(b) denote the result of simulations of the Langevin model with g{k) 
given by 

g{k) = {k-ii)e{k-ii + k2,e)e{ii + k2-k,0), (88) 
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for A = 1, /i = 10, A = 0.941, M = 0.212, A;2 = 4 and ^ = 1. The result shown by 
open squares in Fig. 9(b) deviates from that shown by open circles at | A; — /i |~ 4 where 
effects of multiplicative noises do not work for g{k) given by Eq. (88). It is shown that 
an agreement of Psw{k) with P{k) expressed by open squares is better than that with 
Pqik) expressed by the solid curve (and open circles). 

A similar analysis has been made for the case of p = 0.2, whose results are shown in 
Figs. 10(a) and 10(b). Filled circles and triangles in Fig. 10(a) denote Psw{k) obtained 
by simulations and analytic method, respectively. Solid curve, expressing Pqik) with 
q — 1.20 and a — 1.34, deviates from Psw{k) at | /c — // |~ 4. In contrast, open squares 
express the result calculated with the use of qk given by Eq. (86) for ki — 5 and 9 — 1, 
which are in better agreement with Psw{^^) than Pgik). 

SDE has been apphed to the case of p = 0.2. Open circles in Fig. 10(b) express the 
result of simulations for Langevin equations with g{k) = k — X = 1, fj, = 10, A = 1.788 
and M = 0.111, which are chosen from Eqs. (57) and (59) with q = 1.20 and a = 1.34 for 
p = 0.2 [27]. In contrast, open squares in Fig. 10(b) denote P{k) calculated by simulation 
with g{k) given by Eq. (88) for A = 1, // = 10, A = 1.788, M = 0.111, A;2 = 5 and ^ = 1. 
An agreement of Psw{k) with P{k) is much better than that with Pqik) shown by the 
solid curve (and open circles). 

We have obtained the optimum entropic index of g = 1.0 for p — 0.4, 0.6 and 0.8 
[Fig. 7(c)-7(e)], and q = 0.95 for p = 1.0 [Fig. 7(f)]. Log plots in Figs. 8(c)-8(f) show 
that there are deviations between Psw{k) and Pqik) at | A; — // |~ 5. A Log plot of 
results of p = 1.0 are again shown in Fig. 11, where Psw{k) calculated by simulations 
and analytical method are shown by filled circles and triangles, respectively. We note that 
at I A; — |~ 4, Psw{k) becomes gradually larger than PQ{k) for q — 0.95 and a — 6.0 
shown by solid curve. In order to reduce the discrepancy at | A; — |~ 4, we adopt a 
function oi g{k) for multiplicative noises given by 

^(A;) = e{k-^x-k^,e), (89) 

which has a magnitude of unity a.i k ^ n + k^ and zero cit k ~ ji + k^. This expression 
with Eq. (49) may be alternatively interpreted such that magnitude of additive noises is 
increased by an amount of M at A; ~ // + A;3. Open circles in Fig. 11 show P{k) calculated 
by using SDE with g{k) — k — ji, \ — 1, A — and M = 0, which are chosen from Eqs. 
(57) and (59) with q — 1.0 and a — 2.45 [27]. In contrast, open squares express P{k) 
calculated with g{k) given by Eq. (89) with A = 1, ^ = 6.0, M = 2.0, A;3 = 4 and ^ = 1, 



16 



which is in better agreement with Psw{k) than Pqik) shown by the sohd curve. 

It is necessary to point out that the normahzation factor of l/v^27ro"' in Eq. (69) 
plays an important role in determining 112 (cr') as shown below. When Q-Gaussian of the 
right-hand side in Eq. (69) is replaced by 

P{k) = {k-^)-\ (7>1), (90) 

we get 

/ da e-"^ f(a) = x'^/^ (91) 
Jo 

with l/2cr'^ = a, {k — iif = x, and 

/(«) - ^(D^^Kir). m 

(93) 



V27f \2aJ \\2a 



r(7/2)' 

which lead to 

112(0-) oc a~'^, for large a (94) 

as given by Eqs. (76) and (77). 

On the contrary, if the normalization factor of \j\p2/KG' is neglected in Eq. (76), we 
get Eqs. (91) and (93) with 

/(«) = [^T^M^T\- (95) 



.2a) \\2a 
By using the relation obtained from Eq. (95): 

we get 

n2(cT) oc (7-('^+^\ for large (7 (97) 

The difference between indices of a in Eqs. (94) and (97) arises from the normalization 
factor 1 / \f2n:a' of Gaussian in Eq. (69) . 



7 Conclusions 

In summary, degree distribution of WS small- world networks has been discussed in terms 
of Q-Gaussian obtained by the three approaches: MEM, SDE, and HVD. It is interesting 
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that with an accuracy of | Psw{k) — Pqik) |~ 10^^, a single Pqik) may describe main 
parts of Psw{k) which is expressed as a superposition of Poissons in Eq. (6). In order 
to account for the overall behavior in Ps\v{k) including its main part and tails, we have 
assumed the /c-dependent entropic index in MEM with the entropy given by Eq. (83), 
which is nonextensive unless = 1 for all k. 

This /c-dependent entropic index is realized in SDE, by adopting multiplicative noises 
given by Eq. (88) for p = 0.1 and 0.2, and those given by Eq. (89) for 0.4 ~ p ~ 1.0. 
Degree distributions of scale-free and random networks are obtained for g{k) — k — /i and 
g{k) — 0, respectively, in SDE. We expect that P{k) of SW networks has the combined 
property of scale-free and random networks because g{k) is k — /i and for small \ k — /i \ 
and large | k—ji |, respectively. It is necessary to understand the relation between rewiring 
processes in making WS-SW networks [5] and the two types of multiphcative noises in 
Langevin equation. 

Although the three approaches, MEM, SDE, and HVD, are conceptually rather dif- 
ferent, they lead to equivalent results. There are, however, some differences among them. 

(i) The entropic index q is given by g = 1 + 2M/ (A + M) and g = 1 + 2/72 in SDE and 
HVD, respectively, while that in MEM is a free parameter with the constraint of g < 3. 

(ii) The escort probability P{k) is obtained from as given by Eqs. (22) and (23) in 
MEM, while there are no ways in SDE and HVD. Although no analyses have been made 
for Psw{k) with the use of HVD, Gaussian-based expression (66) or (81) is expected to 
be promising. This subject is left as our future study. 
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Figure 1: Degree distribution Psw{k) (filled marks) of the WS small- world model with 
N = 100 and /i = 10 for various p values, calculated by simulations with 1000 trials, solid 
lines being shown for a guide of the eye. As a comparison, P{k) for random networks 
with the same parameters is plotted by open circles (see text). 

Figure 2: (a) The Q-Gaussian and (b) Pqik) with = 10 and a = 1.0 for various q 
values. 

Figure 3: Time courses of k{t) obtained by dynamical simulations of stochastic differential 
equations [Eqs. (49)-(52)] for (a) A = 1.0 and M = 0.0, and (b) A = 1.0 and M = 0.5 
with A = 1.0 and /v, = 10. Results are simulations of a single trial by using the Heun 
method with the time step of 0.0001. 

Figure 4: (color online) Log-plots of P{k) calculated by simulations (circles) of stochastic 
differential equations [Eqs. (49)-(52)] for (a) A = 1.0 and M = 0.0, and (b) A = 1.0 and 
M = 0.5 with A = 1.0 and = 10, and those of Pq(A;) [Eq. (45)] (solid curves): q = 1.0 
and a = 1.008 for (a) and q = 1.667 and a = 0.984 for(b). Results of simulations are 
averages oi N = 100 for 50 < t < 100 by using the Heun method with the time step of 
0.0001. 

Figure 5: Log plots of 112 (cr) for Q-Gaussian against a/{ao ^pu) with various q values, 
their slopes for g = 1.5, 2.0 and 2.5 being shown by chain curves with 5 = 4, 2 and 1.33, 
respectively (iTo = l/-\/2a^ for g = 1.0). 

Figure 6: The p dependence of R (chain curve), a (dashed curve) and q (solid curve) 
adopted for a fitting of Pqik:) to Psw{k:) of simulations for SW networks with — 100 
and /X = 10, i? — 1 being multiplied by a factor of ten. 

Figure 7: (color online) Pswik) with N = 100 and /i = 10 obtained by simulations 
(circles), and Pqik) (solid curves) with adopted q values for (a) p = 0.1, (b) p = 0.2, (c) 
p — 0.4, (d) p — 0.6, (e) p — 0.8, and (f) p — 1.0. Dashed lines are shown for a guide of 
the eye. 

Figure 8: (color online) Log-plots of Fig. 7. 

Figure 9: (color online) (a) Log plots of Psw{k) for p = 0.1 with fi = 10 and = 100 
calculated by simulations (filled circles) and analytical method (filled triangles), Pqik) 
(solid curve), and P{k) with q^ given by Eq. (86) for ki = 4 and 9 = 1 (open squares), 
(b) P{k) in Langevin method with g{k) = A = 1,A = 0.941 and M = 0.212 (open 

circles), and that with g{k) given by Eq. (88), A = 1, A = 0.941, M = 0.212, A;2 = 4 and 
9 = 1 (open squares) (see text) . 
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Figure 10: (color online) (a) Log plots of Psw{k) for p = 0.2 with n = 10 and = 100 
calculated by simulations (filled circles) and analytical method (filled triangles), -Pq(/c) 
(solid curve), and P{k) with qk given by Eq. (86) for A;i = 5 and 6 = 1 (open triangles), 
(b) P{k) in Langevin method with g{k) = k — n, X = 1,A = 1.788 and M = 0.111 (open 
circles), and that with g{k) given by Eq. (88), A = 1, ^ = 1.788, M = 0.111, A;2 = 5 and 
9=1 (open squares) (see text). 



Figure 11: (color onhne) Log plots of Psw{k) for p = 1.0 with fj, = 10 and N = 100 
obtained by simulations (filled circles) and analytical method (filled triangles), Pqik) 
with q = 0.95 and a = 2.45 (solid curve), P{k) in Langevin method with g{k) = k — fi, 
\ = 1, A = 6.0 and M = and (open circles), and that with g{k) given by Eq. (89), 
\ = 1, A = 6.0, M = 2.0, ks = 4 and 9 = 1 (open squares) (see text). 
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